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Cross-sectional  scanning  tunneling  microscopy  (XSTM)  and  density  functional  theory  have  been  used  to 
characterize  the  structure  of  GaSb  digitally  doped  with  Mn.  The  Mn  dopants  are  found  in  both  isolated 
substitutional  form  as  well  as  in  large  clusters  of  zinc-blende  MnSb  commensurate  with  the  surrounding  GaSb 
matrix.  Theoretical  calculations  predict  that  these  two  forms  of  Mn  in  the  digitally  doped  layers  will  have  a 
very  different  appearance  in  XSTM  images.  Substitutional  Mn  enhances  the  local  density  of  states  near  the 
surface,  thus  appearing  higher  in  a  filled- state  image.  In  contrast,  MnSb  clusters  induce  substantial  structural 
relaxation  at  the  {110}  surface,  and  therefore  appear  as  localized  depressed  regions  with  negligible  perturbation 
of  the  surrounding  GaSb. 
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INTRODUCTION 

Ferromagnetic  materials  that  can  be  easily  integrated  into 
semiconductor  electronics  are  crucial  for  so-called  “spin- 
tronic”  applications.  Concomitant  with  the  obvious  structural 
requirements — for  example,  ferromagnetic  materials  lattice 
matched  to  GaAs  and  Si — is  the  desire  for  room-temperature 
operation.  Diluted  magnetic  semiconductors,  such  as  Mn- 
doped  GaAs,1,2  are  designed  to  fulfill  the  first  of  these  re¬ 
quirements,  whereas  fulfilling  the  requirement  of  room- 
temperature  operation  has  been  more  difficult.  Recently 
some  degree  of  success  in  raising  the  ferromagnetic  ordering 
temperature  has  been  obtained  with  digitally  doped  samples, 
where  the  dopant  is  deposited  periodically  rather  than  con¬ 
tinuously  during  the  growth.3,4  However,  there  has  been  little 
progress  in  understanding  the  atomic-scale  structure  of  the 
doping  profile  in  such  samples. 

In  this  Brief  Report,  we  combine  cross-sectional  scanning 
tunneling  microscopy  (XSTM)  and  density  functional  theory 
to  characterize  the  structure  and  distribution  of  dopants  in 
GaSb/Mn  digital  alloys.  High-resolution  atomic- scale  im¬ 
ages  of  the  (110)  cross-sectional  surface  show  a  number  of 
characteristic  localized  defects.  These  localized  defects  are 
analyzed  using  simulations  of  the  XSTM  experiments  de¬ 
rived  from  density  functional  calculations  and  interpreted  as 
isolated  substitutional  Mn  dopants  on  the  Ga  sublattice,  and 
clusters  of  zinc-blende  MnSb  commensurate  with  the  GaSb 
matrix.  To  the  best  of  our  knowledge,  this  is  the  first  report 
of  a  combined  theoretical-experimental  effort  to  identify 
zinc-blende  MnSb  clusters  in  GaSb. 

EXPERIMENT 

The  sample  was  grown  at  the  University  at  Buffalo  and  its 
magnetoelectronic  properties  have  been  reported  previously.3 
Briefly,  it  was  grown  on  a  p-type  GaAs  (001)  substrate  with 
a  thick,  fully  relaxed  GaSb  buffer  layer.  The  GaSb/Mn  digi¬ 
tal  alloy,  grown  at  275  °C,  consists  of  \  monolayer  (ML)  of 
Mn  deposited  every  16  ML.3  The  sample  was  first  thinned  to 
—200  /xm,  and  then  loaded  into  the  chamber,  scribed,  and 

cleaved  in  situ  to  expose  a  (110)  surface.5  The  XSTM  mea¬ 


surements  were  performed  in  constant-current  mode  in  an 
ultrahigh-vacuum  chamber  with  base  pressure  below  1 
X  10“ 10  Torr.  All  images  shown  here  are  of  filled  electronic 
states,  nominally  showing  the  topography  of  the  Sb  sublat¬ 
tice  on  the  cross-sectional  surface. 

Figure  1(a)  shows  a  nanoscale  overview  of  the  sample, 
including  the  interface  between  the  buffer  layer  and  the  digi¬ 
tally  doped  region;  the  first  18  periods  of  the  GaSb/Mn  digi¬ 
tal  alloy  with  a  period  of  approximately  5  nm  are  indicated. 
Three  different  terrace  levels  are  observed  as  well  as  numer¬ 
ous  steps.  (The  high  step  density  on  the  cleavage  surface  is 
presumably  a  result  of  the  dislocation  network  caused  by  the 
8%  lattice  mismatch  between  the  fully  relaxed  GaSb  buffer 
layer  and  the  GaAs  substrate.)  The  GaSb  buffer  layer  dis¬ 
plays  a  smooth,  uniform  composition;  in  comparison,  filled- 
state  images  of  the  digitally  doped  region  show  frequent 
nanometer-scale  depressions  (dark  areas)  aligned  along  the 
nominal  Mn-containing  layers.  Figure  1(b)  shows  a  higher- 
resolution  image  of  the  region  near  the  interface.  Note  that 
XSTM  of  the  (110)  surface  reveals  only  every  other  row  in 
the  [001]  growth  direction,  so  that  the  observed  period  of  the 
depressions  is  eight  Sb  rows,  consistent  with  the  intended 
Mn  doping  profile  (every  16  growth  layers). 

The  typical  dimensions  of  the  depressions  are  about  1-4 
nm  along  [110]  (in  the  plane  of  the  deposited  layer),  and  1-3 
rows  (0.6-2  nm)  in  the  [001]  growth  direction.  There  is  little 
apparent  perturbation  of  the  lattice  surrounding  these  fea¬ 
tures,  and  no  evidence  of  local  band  bending.  Recent  first- 
principles  calculations  found  that  zinc-blende  MnSb  is  lattice 
matched  to  GaSb  with  an  optimal  lattice  constant  of  6.2  A, 
but  is  mechanically  unstable.6,7  However,  the  theoretical  re¬ 
sults  also  predict  that  a  small  amount  of  Mn  can  be  incorpo¬ 
rated  in  GaSb.6  Considering  that  our  sample  has  been  digi¬ 
tally  doped  with  only  \  monolayer,  we  suggest  that  the 
depressions  observed  in  the  XSTM  images  are  zinc-blende 
MnSb  clusters.  This  is  consistent  with  the  fact  that  the  fea¬ 
tures  appear  to  have  the  same  lattice  structure  as  the  host 
lattice,  rather  than  that  of  pure  Mn  or  MnSb  in  the  NiAs 
structure.  The  0.1  nm  depth  of  the  depressions  could  result 
from  reduced  charge  transfer  from  Mn  to  Sb,  postcleaving 
atomic  relaxation,  or  a  combination  of  such  effects. 
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FIG.  1.  Filled-state  XSTM  images  of  a  GaSb/Mn  digital  alloy 
(2.30  V,  30  pA).  (a)  87  X  87  nm2  area  including  the  end  of  the  GaSb 
buffer  and  the  first  18  superlattice  periods;  (b)  35X35  nm2  area 
indicated  in  (a)  by  the  dashed  square.  Nominal  positions  of  the  Mn 
layers  are  indicated  by  the  arrows.  In  (b)  the  dashed  lines  indicate 
the  region  for  the  line  profile  in  Fig.  5  (see  text). 

In  addition  to  these  localized  depressions  we  observe  a 
number  of  other  characteristic  features;  Fig.  2  shows  an 
atomic-resolution  XSTM  image  of  the  Mn  digitally  doped 
region.  In  addition  to  the  depressions  there  are  several 
atomic-size  bright  and  dark  features  not  normally  observed 
in  GaSb,5,8  suggesting  the  presence  of  additional  foreign  at¬ 
oms  or  defects.  The  small  dark  features,  which  are  present  in 
both  the  GaSb  buffer  layer  and  the  digitally  doped  region, 
are  characteristic  of  As  substitutional  defects9  (substituting 
for  Sb)  and  are  believed  to  be  present  because  the  As  cell 
used  to  prepare  the  GaAs  substrate  was  still  hot  during  the 
GaSb  growth.  As  we  discuss  later  in  the  paper,  the  bright 
features,  which  correlate  strongly  with  the  period  of  the  Mn 
doping,  are  most  likely  due  to  randomly  substituted  Mn  dop¬ 
ants  on  the  Ga  sublattice. 


THEORY 

To  help  interpret  our  experimental  results  we  performed 
electronic  structure  calculations  in  the  generalized  gradient 


FIG.  2.  Filled-state  XSTM  image  of  GaSb  digitally  doped  with 
Mn,  18X18  nm2  (1.90  V,  40  pA).  Arrows  indicate  the  nominal 
positions  of  the  Mn  layers. 

approximation  of  Perdew  and  Wang,10  followed  by  simula¬ 
tions  of  constant-current  XSTM  in  the  approach  of  Tersoff 
and  Hamann.11  A  plane- wave  basis  set  was  used  in  the 
projector-augmented  wave-function  method12  as  imple¬ 
mented  in  the  VASP  code  to  self-consistently  determine  the 
ground-state  charge  density  and  Kohn-Sham  wave  functions, 
with  total  energies  converged  to  0.1  meV.13  The  wave  func¬ 
tions  were  calculated  for  a  single  wave  vector  in  the  Brilloun 
zone  of  the  slab  supercells  at  a  kinetic  energy  cutoff  of  270 
eV.  For  all  the  calculations,  we  used  as  the  lattice  parameter 
our  theoretically  optimized  value  of  6.2  A.  For  the  simula¬ 
tions  of  filled-state  XSTM  images,  the  local  density  of  states 
(LDOS)  was  integrated  over  an  energy  range  of  2  eV  up  to 
the  Fermi  level. 

We  modeled  the  GaSb  (110)  cleavage  surface  in  a  super¬ 
cell  slab  geometry  consisting  of  five  layers  of  GaSb  sepa¬ 
rated  by  a  vacuum  region  of  14  A;  the  bottom  layer  of  the 
slab  was  passivated  to  simulate  semi-infinite  bulk  material. 
We  considered  isolated  Mn  substitutional  defects  in  the  first 
and  second  layers  near  the  surface  and  MnSb  clusters.  For 
the  cases  of  isolated  Mn  dopants  we  used  a  4  X  4  surface  unit 
cell;  for  simulating  clusters  of  MnSb  in  GaSb,  we  used  a  7 
X  1  surface  unit  cell  with  four  bilayers  of  MnSb  in  the  [001] 
growth  direction.  To  minimize  spurious  superlattice  effects, 
for  supercells  containing  a  single  substitutional  Mn,  only  at¬ 
oms  within  7.0  A  of  the  Mn  site  were  allowed  to  relax  during 
the  minimization  of  the  energy.  For  the  supercell  with  four 
MnSb  bilayers,  all  layers  but  the  bottom  passivating  layer 
were  relaxed  during  the  energy  minimization.  The  criterion 
for  convergence  of  the  structural  optimizations  was  that  the 
energy  difference  between  two  structural  configurations  be 
less  than  1  meV. 

Figures  3(a)  and  3(b)  show  the  simulated  filled-state 
XSTM  images  for  Mn  substitutional  dopant  atoms  in  the  first 
and  second  layers,  respectively.  Both  of  these  dopants  have  a 
LDOS  near  Mn  that  is  slightly  larger  than  the  unperturbed 
surface  LDOS,  i.e.,  an  increased  tip-sample  separation.  Spe¬ 
cifically,  these  Mn  substitutionals  increase  the  LDOS  on  sur¬ 
face  Sb  atoms  to  which  they  are  bonded.  Thus  Mn  substitu- 
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FIG.  3.  (Color  online)  Simulated  filled-state  XSTM  images  of 
Mn  substitutional  dopants  near  the  GaSb  (110)  surface.  Ga  and  Sb 
atoms  are  shown  in  gray  (cyan)  and  light  gray  (yellow),  respec¬ 
tively.  Mn  substitutional  dopant  in  (a)  the  first  layer,  and  (b)  the 
second  layer.  A  structural  model  of  the  underlying  lattice  is  shown 
in  (a)  for  the  first  two  layers.  The  location  of  the  Mn  substitutional 
is  indicated  in  each  image. 

tional  in  the  first  layer  is  associated  with  enhanced  LDOS  on 
two  Sb  atoms  [below  the  Mn  site  in  Fig.  3(a)]  and  Mn  sub¬ 
stitutional  in  the  second  layer  is  signified  by  an  enhanced 
LDOS  on  a  single  Sb  surface  atom  [above  the  Mn  site  in  Fig. 
3(b)].  Overall  the  simulated  spectra  are  similar  to  the  bright 
features  observed  in  Fig.  2  and  are  consistent  with  those  of 
substitutional  Mn  dopants  in  GaAs. 14,15  Thus  we  propose  that 
the  atomic-scale  bright  protrusions  observed  in  Fig.  2  are 
associated  with  isolated  Mn  dopants,  and  that  the  larger  hill¬ 
ocks  result  from  overlap  of  the  LDOSs  of  substitutional  Mn 
defects  in  close  proximity. 

To  determine  what  effects  clusters  of  MnSb  commensu¬ 
rate  with  the  GaSb  lattice  would  have,  we  calculated  the 
geometric  and  electronic  structure  of  four  bilayers  of  MnSb 
in  GaSb.  A  structural  model  of  the  optimized  supercell  is 
shown  in  Fig.  4;  note  that  for  clarity,  a  7  X  3  surface  unit  cell 
is  shown.  Although  bulk  MnSb  is  well  lattice  matched  to 
bulk  GaSb,  there  is  substantial  structural  relaxation  near  the 
surface  in  the  MnSb  bilayers.  The  Mn  atoms  at  the  surface 
relax  all  the  way  into  interstitial  locations  below  the  surface 
to  form  complexes  with  Mn  in  the  subsurface  layer.  There  is 
an  accompanying  ~1  A  relaxation  of  the  Sb  surface  atoms 
into  the  surface  in  these  bilayers.  However,  there  is  negli¬ 
gible  relaxation  of  the  surrounding  GaSb.  Thus  small  clusters 
of  MnSb  in  GaSb  should  induce  significant  surface  relax¬ 
ation  on  the  cleaved  {110}  surface  but  little  perturbation  to 
the  surrounding  GaSb. 

The  effects  of  these  structural  relaxations  are  evident  in 
Fig.  5,  where  the  simulated  XSTM  image  of  the  four  MnSb 
bilayers  is  displayed  along  with  a  comparison  of  theoretical 
and  experimental  line  profiles.  The  two  line  profiles  compare 
favorably,  showing  an  approximate  depth  of  1  A,  which 
agrees  well  with  the  predicted  structural  parameters.  In  addi¬ 
tion,  both  the  theoretical  and  experimental  line  profiles  ex¬ 
hibit  similar  structure  in  the  bottom  of  the  depression;  the 
theoretical  results  indicate  that  this  structure  originates  from 
the  surface  Mn  which  relaxes  to  subsurface  locations. 


FIG.  4.  (Color  online)  Top  view  (upper)  and  side  view  (lower) 
structural  model  showing  optimized  coordinates  of  four  bilayers  of 
MnSb  in  GaSb.  Ga,  Sb,  and  Mn  atoms  are  shown  in  gray  (cyan), 
light  gray  (yellow),  and  dark  gray  (red),  respectively.  For  clarity  a 
7X3  surface  unit  cell  is  shown. 


SUMMARY 

Digitally  doped  GaSb/Mn  alloys  have  been  investigated 
with  XSTM  on  the  GaSb  (TlO)  surface  and  interpreted  with 
electronic  structure  calculations.  The  atomic- scale  character¬ 
ization  reveals  coexistence  of  isolated  Mn  substitutional  dop¬ 
ants  on  the  Ga  sublattice  and  quasi-two-dimensional  islands 
of  zinc-blende  MnSb  in  the  digitally  doped  layers.  Atomic 
relaxation  of  the  (110)  cleavage  surface  on  the  MnSb  islands 
produces  Mn  substitutional-interstitial  complexes  with  dis¬ 
tinct  signatures  in  the  XSTM  images. 


FIG.  5.  (Color  online)  Simulated  filled-state  XSTM  results  for 
four  bilayers  of  MnSb  in  GaSb.  Ga,  Sb,  and  Mn  atoms  are  shown  in 
gray  (cyan),  light  gray  (yellow),  and  dark  gray  (red),  respectively, 
(a)  Simulated  STM  image  for  the  (110)  surface  with  the  top  few 
atomic  layers  projected  onto  the  data  for  clarity;  (b)  comparison  of 
experimental  and  theoretical  line  profiles  through  the  MnSb  bilay¬ 
ers.  In  (b)  the  zero  of  the  vertical  scale  is  arbitrary.  The  experimen¬ 
tal  line  profile  is  from  the  dashed  region  indicated  in  Fig.  1(b). 
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